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ABSTRACT 

We present the results from 10 nights of observations of the globular cluster NGC 698 1 (M72) 
in the V, R and / Johnson wavebands. We employed the technique of difference image anal- 
ysis to perform precision differential photometry on the time-series images, which enabled us 
to carry out a census of the under-studied variable star population of the cluster. We show that 
20 suspected variables in the literature are actually non-variable, and we confirm the variable 
nature of another 29 variables while refining their ephemerides. We also detect 1 1 new RR 
Lyrae variables and 3 new SX Phe variables, bringing the total confirmed variable star count 
in NGC 6981 to 43. 

We performed Fourier decomposition of the light curves for a subset of RR Lyrae stars 
and used the Fourier parameters to estimate the fundamental physical parameters of the stars 
using relations available in the literature. Mean values of these physical parameters have al- 
lowed us to estimate the physic al parameters o f the parent cluster We derive a metallicity of 
rFe/H17^ y ~ -1.48±0 03 on the IZjniL&Wed (Il9841) scale (or [Fe/H]uvEs « -1.38±0.03 on 
the new lCarretta et all ( l2009l) scale) for NGC 6981, and distances of ~16.73±0.36 kpc and 
~16.68±0.36 kpc from analysis of the RRO and RRl stars separately. We also confirm the 
Oosterhoff type I classification for the cluster, and show that our colour-magnitude data is 
consistent with the age of ^12.75±0.75 Gyr derived by Pott er et al.. (2010). 
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1 INTRODUCTION 

The study of Galactic globular clusters is important for many rea- 
sons. These stellar systems represent some of the oldest, and con- 
sequently metal poor, stellar populations in the Galaxy, and their 
scrutiny allows us to glean information regarding the formation and 
early evolution of the Galaxy. The spatial distribution of the clusters 
reveals a different aspect of the Galactic structure than other stars in 
the Galaxy, and their orbits and tidal tails provide constraints on the 
Galactic potential. Of course, what we learn about globular clusters 
in our own Galaxy is also applicable to other galaxies as well. 

Globular clusters are also believed to be a close approximation 
to a stellar laboratory since a cluster's members were formed at the 
same time from the same primordial material with the same compo- 
sition, leading to a homogeneity of certain properties within each 
cluster, but with differences in these properties between clusters. 
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Although this paradigm is being challenged by the recent discovery 
in some globular clusters of multimodal main sequences and sub- 
giant branches (Piotto 2009, and references therein), indicating the 
existence of multiple stellar populations, most globular clusters do 
not exhibit such obvious deviations from a simple stellar population 
and the paradigm still holds. 

There are ~ 1 50 globular clusters in our Galaxy for which their 
fundamental properties, such as metallicity, dist ance, age and kine- 
matics, have been estimated by various methods faarrisll993h . One 
independent method for estimating at least some of these quantities 
is by studying the population of RR Lyrae variable stars present 
in most clusters. This method uses the fact that the light curve 
morphology of RR Lyrae stars is connected with their fundamental 
stellar parameters, and consequently quantities such as metallicity, 
absolute magnitude and effective temperature may be calculated 
from the fit parameters of the Fourier decomposition of their light 
curves using empirical, semi-empirical or theoretical relations pub- 
lished in re cent years jSimon & Clemen ^1 19931: lJurcsik & Kovacsl 
19961: Ijures ik 199 8: iKovacj 1 19981 : iKovacs & w5ke3 1200 ll : 
Morgan. Wahl & Wieckhorstll2007h . Appropriate mean values of 
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these fundamental parameters then enable similar estimates of the 
physical parameters of the parent cluster. 

As part of a series of papers on detecting and characterising 
the variable stars in globular clusters, and usi ng the results to esti- 
mate th e parameters o f the hos t cluster ( Arellano Ferro et alj 



2004 iLazaroetalJ |2006|; [Arellano Ferro et alj 120081: 



Arellano Ferro et ai] l2008bl : [Arellano Ferro et ai] l2010h 



have performed CCD time-series photometry of the globular 
cluster NGC 6981 (RA a = 20''53"'27.9'\ Dec. 5 = -12°32'13", 
J2000; / = 35.16°, b = -32.68°) using the method of difference 
image analysis (Section The known variables in this cluster, 
which are exclusively RR Lyrae variables, hav e been studied in a 
handful of photographic observing campaigns ( Shapley & Ritchie 



handtul or photographic observmg campaigns ( shapley & Kitchie 
1923; lRosindll953l : ISawveilll953l : lNobililfl957. : Dickens & Flinr 
1973), the most recent of which is now 40 years in the past. 



Periods and ephemerides are poorly determined or non-existent 
for a substantial number of variables, and light curves for many 
of the claimed "variables" have not been published. Therefore, in 
Section[3l we use our precision differential photometry to perform 
an essential variable star census for the cluster. 

In Section|4l we use the Fourier decompostion of the RR Lyrae 
star light curves to estimate their fundamental physical parameters, 
and then in Section |5] we use the RR Lyrae properties to estimate 
the metallicity of, and distance to, NGC 6981. We also discuss the 
age estimates for the cluster that are available in the literature in 
the context of our colour-magnitude diagram. Our conclusions are 
presented in Section|6] 

Through out this paper we adopt the RR Lyrae nomenclature 
introduced by Alcock et al.l ^2000): namely, RRO designates an RR 
Lyrae star pulsating in the fundamental mode, and RRl designates 
an RR Lyrae star pulsating in the first-overtone mode. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Observations 

We employed the 2.0m telescope of the Indian Astronomical Ob- 
servatory (lAO), Hanle, India, located at 4500m above sea level, 
to obtain time-series imaging of the globular cluster NGC 6981. 
The image data were obtained during several runs between October 
2004 and September 2009, where we collected a total of 103, 110 
and 3 images through Johnson V , R, I filters, respectively (see Ta- 
ble[T]for a detailed log of the observations). The CCD camera that 
was used is equipped with a Thompson CCD of 2048x2048 pix- 
els with a pixel scale of 0.296 arcsec pix~^ and a field-of-view of 
-10.1x10.1 arcmin^ 



2.2 Difference Image Analysis 

As in previous papers i Arellano Ferro et al. 2004 \_ Lazaro et al.l 



2006'; 'Arellano Ferro et al J l2008al ; [Arellano Ferro et all l2008bl : 



Arellano Ferro et al. 2010j), we have employed the technique of dif- 
ference image analysis (DIA) to extract high precision photometry 
for all point sources in the images o f NGC 698 1 , in cluding those 
in the highly crowded cen tral region l lAlard & Lupto n 1998; Alard 
I2OOOI ; iBramich et al. 2005). We used a pre-release version of the 
D a n D 1 43 pipeline for the data reduction process (Bramich et al., in 



^ DanDIA is built from the DanlDL library of IDL routines available at 
|http://www.danidl.co.uk| 



Table 1. The distribution of observations of NGC 698 1 for each filter, where 
the columns Ny , Nn and Nj represent the number of images taken for the 
filters V, R and /, respectively. We also provide the exposure time, or range 
of exposure times, employed during each night for each filter in the columns 

tY,tji and tj . 



Date 


Nv 


tv (s) 


Nr 


tR (S) 


N, 


ti (s) 


20041004 


20 


60-180 


19 


100-150 







20041005 


29 


90 


29 


70 







20050514 


12 


150 


15 


120 







20050515 


13 


150 


15 


120 







20050516 


5 


100-150 


6 


80-120 







20070522 


4 


200-300 


5 


200-240 







20070804 


5 


240 


8 


180 







20070905 


11 


150 


11 


120 







20090913 


4 


150 


1 


120 


3 


100 


20090914 







1 


120 







Total: 


103 




110 
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preparation) which includes a new algorithm that models the con- 
volution kernel matching the point-spread functio n (PSF) of a pai r 
of images of the same field as a discrete pixel array teramicj2008h . 

In brief, the DanD I A pipeline was used to perform the follow- 
ing processing steps on the raw image data. Bias level and flat field 
corrections were applied to the raw images to create calibrated im- 
ages, which were then cleaned of cosmic r ays using the cosmic ray 
cleaning algorithm of lvan Dokkunj JioOlh . A reference image for 
each filter was constructed by registering and stacking a set of the 
best-seeing calibrated images such that all images used were taken 
on a single night. This resulted in 7, 7 and 2 images being stacked 
with total exposure times of 630, 490 and 200 s for the filters V, R 
and /, respectively. The full-width half-maximum (FWHM) of the 
PSF in the V, R and / filter reference images was measured to be 
~3.8, ^--^3.6 and ~4.7 pix, respectively. 

In each reference image, we measured the fluxes (referred to 
as reference fluxes) and positions of all PSF-like objects (stars) by 
extracting a spatially variable (with polynomial degree 3) empirical 
PSF from the image and fitting this PSF to each detected object. De- 
blending of very close objects was attempted. Using the Delauna y 
triangulation and triangle matching method of lPal & BakosI ( l2006h . 
the detected stars in each image in the time-series sequence were 
matched with those detected in the corresponding reference image, 
and a linear transformation was derived which was used to register 
each i mage with the r eference image using cubic O-MOMS resam- 
pUng ( IbIu et al.ll200ll) . 

For each filter, a sequence of difference images was created 
by subtracting the relevant reference image, convolved with an ap- 
propriate spatially variable kernel, from each registered image. The 
spatially variable convolution kernel for each registered image was 
determined using bilinear interpolation of a set of kernels, mod- 
elled as pixel arrays, that were derived for a uniform 6x6 grid of 
subregions across the image. 

The differential fluxes for each star detected in the reference 
image were measured on each difference image as follows. The 
empirical PSF at the measured position of the star on the reference 
image was determined by shifting the empirical PSF model cor- 
responding to the nearest pixel by the appropriate sub-pixel shift 
using cubic O-MOMS resampling. The empirical PSF was then 
convolved with the kernel model corresponding to the star posi- 
tion and current difference image, and then optimally scaled to the 
difference image at the star position using pixel variances a^^j for 
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Figure 1. Plot of standard magnitude in the Johnson system against mean 
instrumental magnitude for the set of standard stars supplied by P. Stetson 
that lie in our field of observation (open circles). The top and bottom pan- 
els coiTespond to the V and I filters, respectively. The solid fine in each 
panel shows the fitted relation between the standard and mean instrumental 
magnitudes. 



image k, pixel column i and pixel row j, taken from the following 
standard CCD noise model: 



2 



(1) 



where ctq is the CCD readout noise (ADU), Fij is the master flat- 
field imag e, G is the CCD gain (e~/ADU) and Mkij is the image 
model (see lBramictJl2008l) . 

Light curves for each star were constructed by calculating the 
total flux /lot (t) in ADU/s at each epoch t from: 



/tot(t) = /rcf + 



/d.ff(t) 
Pit) 



(2) 



where /rcf is the reference flux (ADU/s), /diff(t) is the differential 
flux (ADU/s) and p(t) is the photometric scale factor (the integral 
of the kernel solution). Conversion to instrumental magnitudes was 
achieved using: 



m,„,(t) = 25.0 -2.5 log(/,ot(t)) 



(3) 



where mins(t) is the instrumental magnitude of the star at time t. 
Uncertainties were propagated in the correct analytical fashion. 



2.3 Caveats Of Difference Imaging Analysis 

The value of the reference flux /,cf for each star is measured on the 
reference image, and where the star field is very crowded, the mea- 
sured /ref values are likely to be systematically too large due to flux 
conta mination from blending and unmo delled faint background 
stars jTodd et al.|[2005l: IXodd et al.ir200d) . On the other hand, the 
measurement of the difference fluxes /diff(t) does not suffer from 
this problem because the majority of sources are fully subtracted in 
the difference images, including the blended stars and faint back- 
ground stars. For a variable object, a value for /,cf that is system- 
atically too large will result in a light curve with an amplitude in 



magnitudes that is systematically too small, and viceversa (Equa- 
tions |2]&[3](, although the "shape" of the light curve will be unaf- 
fected. Consequently, we have made a note in column 4 of Table |2] 
(see later) of those variable stars likely to be affected by flux con- 
tamination, or that are affected by a blended PSF. 

Furthermore, it is worth mentioning a feature of the difference 
image construction with DanDIA which can impede photometry 
on some objects. Saturated pixels in the reference image are flagged 
as bad pixels, and when the reference image is convolved with a 
kernel in order to match the PSF of the current registered image, 
the bad pixels in the reference image are grown by the footprint of 
the convolution kernel. For our data, we modelled the kernel as a 
circular pixel array of radius equal to twice the PSF FWHM in the 
current registered image. Hence a saturated star in the reference im- 
age discounts an area in the difference images that encompasses the 
saturated star and its immediate neighbourhood, with poorer see- 
ing resulting in a larger discounted area. Consequently stars in the 
neighbourhood of saturated stars may suffer from imprecise pho- 
tometry (because fewer good pixels are available in the difference 
images for photometric measurements), missing photometric mea- 
surements for a subset of epochs (when the seeing is too poor), 
or simply it may have been impossible to extract any photometric 
measurements at all for the star in question. 

We made our choice of images to be combined into the stacked 
reference images so as to minimise the number of saturated stars. 
However, our reference images still contain a handful of saturated 
stars towards the centre of the cluster, which has affected the pho- 
tometry of some of the variable stars. We will refer to this point a 
number of times later in this paper. 



2.4 Photometric Calibrations 

We derived photometric calibration relations for the conversion of 
instrumental v a nd i magnitude s to the Johnson-Kron-Cousins pho- 
tometric system ( iLandoltl l 19 92b by analysing the standard V and / 
magnitudes of a set of standard stars in the observed field of the 
cluster which were kindly made available to us by P. Stetson (pri- 
vate communication). We identified 26 and 28 standard stars in the 
V and I reference images, respectively, and fitted a linear relation 
between standard and mean instrumental magnitudes for each fil- 
ter (see Figure[TJ. No colour term was found to be significant, and 
the linear correlation coefficients are >0.999. These relations were 
then used to convert the instrumental v and i photometric measure- 
ments for all detected point sources to the standard Johnson-Kron- 
Cousins photometric system. 

In Figure |2(a)| we plot the root-mean-square (RMS) magni- 
tude deviation for each of the 8199 calibrated V light curves versus 
the mean magnitude. We achieve better than 20 mmag scatter for 
stars in the magnitude range 14 to ~18.5 mag, and ~6-10 mmag 
photometry at the bright end. 

The instrumental r magnitudes were retained in the instru- 
mental system since no standards with R magnitude measurements 
were found in the literature. In Figure [2(b)] we plot the RMS mag- 
nitude deviation for each of the 9260 instrumental r light curves 
versus the mean magnitude. In this filter we achieve ~8-12 mmag 
photometry at the bright end. 



2.5 Astrometry 

A linear astrometric solution was derived for the V filter reference 
image by matching ~1000 hand-picked stars with the USNO-Bl.O 
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Figure 2. (a) Plot of RMS magnitude deviation versus mean magnitude for each of the 8199 cahbrated V Hght curves, (b) Plot of RMS magnitude deviation 
versus mean magnitude for each of the 9260 instrumental r light curves. Both plots: RR Lyrae variables of the RRO and RRl types are plotted as solid and 
open blue circles, respectively. Variables of the SX Phe type are plotted as solid red circles, and suspected RR Lyrae variables are plotted as solid green circles. 
Previously suspected variables in the literature that do not show any variability in our data are marked as red crosses. A few of these previously suspected 
variables exhibit a relatively large RMS magnitude deviation which is due to outher photometric measurements in the light curve rather than bona-fide 
variability. 



star catalogue jMonet et a l."2003') using a field overlay in the im- 
age display tool GAIA ([Draper 2000). We achieved a radial RMS 
scatter in the residuals of ^^0.3 arcsec. The astrometric fit was then 
used to calculate the J2000.0 celestial coordinates for all the con- 
firmed variables in our field of view (see Table |3]l. The coordinates 
correspond to the epoch of the V reference image, which is the 
heliocentric Julian date ^2453284. 11 d. 



3 VARIABLE STARS IN NGC 6981 

The first clai med detectio n of two variable stars in NGC 6981 
was made bv lDavisI ( Il917h . However, the first proper study of the 
variab le stars in NGC 6981 was undertaken bv IShaplev & Ritchie! 
l ll920l) (from now on SHI 920) where 18 photographic plates taken 
over 2 years from the 60-inch reflector at the Mount Wilson Ob- 
servatory were analysed. Stars V1-V34 were listed as "probably 
variable" by Miss Ritchie and photometric measurements were ex- 
tracted relative to a set of 29 comparison stars of "sensibly constant 
light". However, it was noted that stars V6, V19, V22, V26 and 
V33 "do not appear to be conspicuously variable", and that V25, 
V30 and V34 "undoubtedly vary, but it has not been possible to 
obtain uniform periods for them". Hence SH1920 actually claim 
detectio n of 29 vari ables, and provide periods for 26 of them. 

Rosind studied 22 of the claimed variables from the 

work of SH19 20 and estimate d new periods for 16 of them. Around 
the same time. lSawveil ( ll953h claimed the detection of 7 more vari- 
ables (stars V35-V41) from a collection of 61 photographic plates, 
although no light curves were published and no periods were es- 
timated. [Nobili ( 1957} further studied the stars V7, V29 and V35, 
presenting lig ht curves and new period s. 

Later on. lDickens & FlinrJ i ll 9721) analysed B and V filter ob- 
servational data for 21 RR Lyrae stars with clear variations, using 
~20 photographic plates for each filter taken on 13 nights spread 
over the period of ~1 year. They revise the periods of these stars 
by combination with previous data, and attempt to use the RR 
Lyraes to derive variou s cluster properties. From the same data, 
iDickens & Flinnl ( Il973h claim the discovery of "a new red vari- 
able", but no light curve or period was ever published. This star 



was subsequently labelled as V42 in the lSawver Hoed i ll 9731) cat- 
alogue. 

In a brief paper bv lKadla et al.l ( Il995l) . the authors use the clus- 
ter colour magnitude diagram to identify 16 stars in the RR Lyrae 
instability strip that are not already claimed to be variable in pre- 
vious publications. The stars, labelled S1-S9 and R1-R7, are put 
forwards as "suspected" and "possible" variables, respectively. Al- 
though these stars are not confirmed variables, they would be good 
candidates to look for variability if time-series photometric data 
exists. 

As we have mentioned in the introduction, no time-series pho- 
tometry for NGC 6981 has been published for the last ^^40 years, 
and none of the time-series photometry that has been published was 
obtained using modern CCD imaging cameras. Our time-series ob- 
servations provide more data points (~100 compared to ~20-60) 
over a longer time-base (~5 yr compared to ~l-2 yr) than the pho- 
tographic campaigns of previous authors, which gives us the po- 
tential to detect variability on timescales from hours to months. 
We also reach deeper by ~4 magnitudes, and achieve a much 
better photometric precision per data point ( ^ 20 mmag down to 
~18.5 mag compared to ~50 mmag at ~17 mag), which is mainly 
the consequence of the use of a larger telescope coupled with a sen- 
sitive CCD imager. Also, the technique of difference imaging has 
enabled us to detect variables in the most crowded central parts of 
the globular cluster. Therefore, we are in a position to fully revise 
the list of variables in the cluster. 

In Table [2] we present the details of all the confirmed vari- 
ables in NGC 6981; namely, those variables for which we have 
detected brightness variations in our light curves above the noise 
level, or which have clearly variable light curves published in the 
literature. The latter of these criteria covers the cases of the two 
RR Lyrae variables V27 and V35 that lie outside of our field of 
view, and which therefore do not have a light curve in our data (see 
Section [ni l. We also provide a comprehensive set of celestial co- 
ordinates for all confirmed variables in our field of view in Table|3] 
and a set of finding charts in Figure|3] since clear finding charts dis- 
playing all confirmed variables have not previously been published. 
This publication should therefore serve as the definitive reference 
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Table 2. Details of all confirmed varia bles in NG C 698 1 . The variable star "V" identification numbers are listed in the first column, and the previous star 
identification numbers (for the Kadla e t alj candidate variables) are listed in the second column. The third column indicates the variable star type. A 

guide to the level of blending of the variable star PSF in the reference images is provided in column 4, where "Inner" and "Outer" indicate that the vaiiable 
is inside or outside, respectively, of the highly crowded central area of the cluster (the dividing line is placed at r = 50 arcsec), and where "Blend" indicates 
that the vaiiable is blended with a star of similar or greater brightness. In column 5, we list the number A^max and A^min of light curve maxima and minima, 
respectively, that were observed in full during our observation runs. Note that the epoch of maximum light T,„ax in column 6 is a heliocentric Julian date (HJD). 
To enable easy referencing, the period for each variable as found by previous authors is supplied in columns 8 to 11 , and may be compared to our derived 
periods in column 7. Note that most of the period estimates from previous authors fail to phase our light curves properly, which highlights the importance of 
the longer temporal baseline and much better photometiic precision of our data. 
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""These variables are RR Lyrae stars, and we believe that they are of the RRO type, but we cannot confirm this (see Section [X6l . ''The classification for this 
vaiiable is most likely RRO (see Section[33). "^These variables are RR Lyrae stars, but we have been unable to distinguish their sub-type (see Section [X6l . 
'*The epoch of maximum light is uncertain because we have not observed the light curve peak. The epoch reported here is that of the data point closest to the 
suspected peak. "^The period listed for V29 is the period Pq in Equation|3] We detect a secular period change for this star of /3 — 1.38 X 10^** d d^^ (see 
Section[33)- -^Due to the poor phase coverage of our observations for this variable, we have been unable to determine a reliable period (see Section lXS^ . 
9The periods listed for the SX Phe variables are the periods associated with the largest amplitude oscillation. '"Our photometric data are not well-phased by 
these periods. 
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Figure 3. Finding charts constructed from our V reference image; north is up and east is to the right. The cluster image is 5.31' by 4.93', and the image stamps 
are of size 23.7" by 23.7". Each confirmed variable lies at the centre of its corresponding image stamp and is marked by a cross-hair. Note that finding charts 
for V27 and V 35 are not available from our data because these variables are not in our field of view. The best finding chart for these stars may be found in 
iDickensI Ii9t3> . 
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Table 3. Celestial coordinates for all the confirmed variables in our field of view, except V27 and V35 which lie outside of our field of view. The coordinates 
correspond to the epoch of the V reference image, which is the heliocentric Julian date ^2453284. 11 d. 
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Table 4. Time-series V , r and / photometry for all the confirmed variables in our field of view, except V27 and V35 which lie outside of our field of view. 
The standard Mstd and instrumental mi„s magnitudes are listed in columns 4 and 5, respectively, coiTesponding to the variable star, filter, and epoch of mid- 
exposure listed in columns 1-3, respectively (note that we do not supply standard R magnitudes). The uncertainty on m\^^ is listed in column 6, which also 
corresponds to the uncertainty on Afsij. For completeness, we also list the quantities /,ef, /jiff and p from Equation |2] in columns 7, 9 and 11, along with 
the uncertainties a^^f and crjiff in columns 8 and 10. This is an extract from the full table, which is available with the electronic version of the article (see 
Suppoiling Information). 
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for the known variable star population of NGC 698 1 , and the reader 
should not need to refer elsewhere in order to collate the pertinent 
information on these variables. 

Our V , r and / time-series photometry for all the confirmed 
variables in NGC 6981 in our field of view is available in Table|4l 
of which we reproduce only a small part in this paper, and which is 
available in full in electronic form (see Supporting Information). 

The following sections describe our findings for the variables 
in NGC 6981, including our detection of previously unknown vari- 
ables. 

3.1 Stars Without Light Curves 

The stars V27, V35, V38 and V42 do not have light curves in our 
data set. The reason for this is that V27 and V35 are not in our field 
of view, and V42 is saturated in our reference image for each fil- 
ter. For V38 we were unable to perform photometric measurements 
on the difference images because the star lies very close to a satu- 
rated star in our reference image for each filter (see Section |273l l. 
We note that both V27 and V35 have light curves published in 



Dickens & Flinnj ( Il972l) . and V35 has a light curve published in 
iNobilil (il957i) . These published light curves show clear variations 
of the RRO type, and periods for these variables have been mea- 
sured in more than one publication. Consequently we include these 
stars as confirmed variables in Table|2] 

The stars V38 and V42 have not been included in Table |2] 
because there have been no published light curves or periods for 
these two stars, and the only claim to their v ariab ility are the brief 
statem ents of this property in I Saw veil ( Il953l) and lOickens & Flim] 
l ll973h . The proximity of V38 to a much brighter star m ay have 
confused its identification as a variable bv lSawvej ( [l953l) . There- 
fore, further time-series observations will be required to promote 
either of these stars to confirmed variables. 

3.2 Stars That Do Not Show Variability 

Considering first the stars V1-V34 fisted in SHI 920, we find that 
our fight curves for V6, V19, V22, V26, V30, V33 and V34 do 
not show any variability, and that the ~7-20 mmag RMS scatter 
in the light curves is consistent with the noise for the majority of 
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constant stars of '--^ITth magnitude (see Figure|2ll. This confirms the 
comment made in SH1920 that V6, V19, V22, V26 and V33 "do 
not appear to be conspicuously variable", and refutes the statement 
that V30 and V34 "undoubtedly vary". In the light of our precision 
photometry we have dropped these stars from the variables listed 
in Table|2] 

Similarly, we find that V37, V39, V40 and V41 , from the list 
of stars V35-V41 claimed to be variable bv lSawveJ l fT953h . show 
no variability in our light curves, although the light curves of V37 
and V39 suffer from some outlier photometric measurements lead- 
ing to artificially large RMS magnitude deviations (see Figure O. 
Hence, we have also dropped these stars from the variables listed 
in Table|2] 

Finally , for th e variable star candidates S1-S9 and R1-R7 from 
iKadla et al 1( ll995l) . we do not detect any variability above the ~8- 
35 mmag light curve noise level for the stars SI, S2, S4, S6, Rl, R2, 
R5, R6 and R7, and therefore we have dropped these stars from the 
variables listed in Table [2] The remaining variable star candidates 
are RR Lyrae stars (see Section [33T l. 

3.3 Detection Of New Variable Stars 

Stars V43-V56 are newly detected variables. Of these, V43-V49 
were discovered by inspection of the light curves of the variable star 
candidates from lKadla et al.l il995h . and V50-V53 were identified 
in the sequence of difference images as clearly varying sources. 
By converting each difference image Duij to an image of absolute 
deviations in units of sigma -D^ij = |-Dfeij l/f^ij (see Equation[TJ 
and then constructing the sum of all such images Sij = ^fcij 
for each filter, one can identify candidate variable sources as PSF- 
like peaks in the image Sij. Using this method we discovered the 
new variables V54-V56 and we recovered all of the other variable 
stars listed in Table [2] 

We also searched for variability in the set of light curves by 
apply ing the string-length method (Burke etal. 197G; Dworetsky 
Il983l) to each light curve to determine the best period and a corre- 
sponding normalised string-length statistic Sq. Exploration of the 
light curves with the smallest values of Sq by phasing them with 
their best period recovered most of the RR Lyrae stars in our list 
of variables, but did not yield any further variable star detections. 
In conjunction with the string-length search method, we inspected 
the light curves of all stars with an unusually large RMS magnitude 
deviation using the plots in Figure[2]as a guide to the typical RMS 
values at a given magnitude. Again, we did not discover any more 
variable stars. 

A useful result of the string-length period search was that we 
had obtained initial period estimates for the variable stars, from 
which we were able to determine accurate periods by performing 
another string-length period search around the initial period esti- 
mates, but using a finer grid of test periods. Care was taken with 
the analysis of light curves showing the Blazhko effect (see Sec- 
tion [331( and in these cases we used appropriate subsets of self- 
consistent data to obtain the correct period. Note that we were un- 
able to determine periods for the variables V44, V45, V52 and V53, 
and that the period for V51 is unreliable (see Section [33] ). For the 
SX Phoenicis stars V54-V56, we used a different method to deter- 
mine their periods (see Section lSTTt . 

3.4 Variable Stars With Blended PSFs 

NGC 6981 has a very crowded central concentration of stars which 
implies that many of the variables we have observed are blended 



with other stars, and those that we observed towards the centre of 
the globular cluster are almost certainly blended with one or more 
sources. In most cases the blend stars are much fainter than the 
variable stars (which are dominated by bright RR Lyrae stars) and 
will have little effect on the measured reference fluxes and light 
curve variation amplitudes (see Section |23] ). However, as a useful 
guide to the likelihood of a variable including a blend, in column 4 
of Table [2] we indicate whether a variable star lies in the crowded 
centre of the cluster or not, by setting a dividing line at a radius of 
50 arcsec within which no sky background is visible in our refer- 
ence images. 

The PSF of all confirmed variable stars in the reference images 
was inspected in detail in order to identify blended PSFs. We find 
that V18, V24, V36, V45, V51 and V52 are closely blended with 
stars of similar or greater brightness, and in these cases the vari- 
able star reference flux may be unreliable, depending on the suc- 
cess of the deblending iterations applied by the DanDIA software. 
The blend status of these variables is clearly marked in Table[2l 

3.5 RR Lyrae Stars 

In Table [2] we present data on 36 confirmed RR Lyrae stars, 31 
of which we classify as RRO type, and 5 as RRl type, although 
the classification for V51 is a little uncertain. Classification was 
based primarily on light curve shape and period, and confirma- 
tion of our results was achieved by constructing the Bailey diagram 
of amplitude versus period, which clearly separates the RRO stars 
from the RRl stars (see Figure|6l(. In the Bailey diagram, we also 
plot the ave rage distribution of the RRO and RRl variables in M3 
dCacciari. C orwin & Carnev 2005), the prototype Oosterhoff type I 
cluster, for comparison purposes. We note that V27 has an anoma- 
lously large amplitude (taken from lPickens & Flinnlll973) for its 
period, and that V46 has an anomalously small amplitude for its 
period. 

Phased V and r light curves of the 29 RRO and 5 RRl stars for 
which we have data are displayed in Figures|4]&|5] respectively. Of 
these variables, 5 RRO and 2 RRl are new discoveries (V47-V51 
and V43,V46). 

Our observations for V5 1 have poor phase coverage such that 
they only cover the falling part of the light curve. Consequently, 
we have been unable to derive a reliable period, and the period 
P = 0.357335 d that we present in Table [2] derived using the 
string-length method is probably a smaller period alias of the real 
period. There is no information in the light curve that enables us 
to rule out other plausible period aliases. Although V51 is blended 
in the reference image, it still lies neatly in the RRO dominated re- 
gion of the instability strip in the horizontal branch of the V — r 
colour-magnitude diagram (CMD; see Figure|7}. Based on this evi- 
dence, V5 1 is clearly an RR Lyrae star, most likely of the RRO type. 
Further evidence for the RRO type comes from the non-sinusoidal 

shape of the light curve at the derived period. 

^It is currently believed that the Blazhko effect telazhkd 

Il907h . a periodic modulation of both the light curve amplitude 
and phase on timescales of typically tens to hundreds of days, oc- 
curs frequently ( > 40 per cent) in RRO stars tJurcsik et al.l l2009l : 
iKolenberg etal] l2010h . Until recently the incidence rate of the 
Blazhko effect in RRO stars was estimated to be substantially lower 
(e.g. ~25 per cent in Smith 1981). However, this picture has be- 
come clearer with the execution and analysis of long-time-base 
quasi-continuous photometric observations of high precision en- 
abling one to detect smaller amplitude manifestations of the effect, 
and to detect the cases of RRO stars exhibiting unstable Blazhko 
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Figure 4. Phased V and r light curves of the 29 RRO stars for which we have data, using the periods listed in Table|2] The light curve for V29, marked with an 
asterisk, is phased with the period and secular period change detailed in Table|2] For V32, the only data to show the Blazhko effect are the observations from 
the nights in 2005, which do not phase well with the rest of the light curve, and are marked in the corresponding plot panel by vertical-dashed data points. 
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Figure 5. Phased V and r light curves of the 5 RRl stars for which we have data, using the periods listed in Table[2] 



Bailey Diagram : NGC 6981 



0,5 




V27 



-0.6 



og 



Figure 6. Bailey diagram of the light curve V amplitude versus (log-)period 
for the RR Lyrae stars in NGC 6981. Variables of the RRO type are plotted 
as solid circles where the light curve amplitude is well-defined between 
observed maxima and minima, and they are plotted as open triangles oth- 
erwise (which then represent a lower bound to the Y amplitude). Variables 
of the RRl type are plotted as open circles. V51 has been omitted from 
this plot because it does not have a reliable period estimate. The continuous 
lines represent the average distribution o f the RRO and RRl variables in M3 
from lCacciari, Corwin & Carney] j2005h . The clear outliers in this diagram 
are labelled with their variable star designations, namely V27 and V46. 
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Figure 7. A portion of the V — r CMD for NGC 6981, focussing on the 
horizontal branch and part of the red giant branch. RR Lyrae variables of the 
RRO and RRl type are plotted as solid and open blue circles, respectively. 
Three of the suspected RR Lyrae variables are plotted as solid green circles. 



effects such as abrupt changes in the modulation ampUtude and/or 
the appearance/disappearance of the Blazkho effect. Note that the 
Blaz hko effect has also been dete cted in a number of RRl stars 
(e.g. lWils, Kleidis & Broeiisll2008h . 

Our observations cover a ~5 yr baseline with intense obser- 
vations during small groups of nights. The ^^10 mmag precision 
that we have achieved at 1/ =17 mag allows us to detect small 
changes in amplitude and/or period of the RR Lyrae stars on the 
typical timescale of the Blazhko effect. The RRO variables VII, 
V23, V28, V31 and V32 show a Blazhko effect that is clearly visi- 
ble in the phased light curves in Figure|4] Further careful inspection 
of the phased light curves (at a smaller scale and using different 
plot symbol colours for each observation night) has allowed us to 
detect a small amplitude Blazhko effect in the light curves of the 



RRO variables VIO, V14, V15, V36 and V49, which is not so visi- 
ble in the phased light curves in Figure|4]due to the large scale that 
has been used. A more detailed analysis of the periodicity of the 
Blazhko effect in these variables is not possible with our data since 
we have observed very few light curve maxima (see column 5 of 
Table I2J- Note that we do not find evidence for the Blazhko effect 
in the light curves of the RRl variables. 

Having detected the Blazhko effect in 10 out of 29 RRO vari- 
ables, we can safely place a ~34 per cent lower limit on the rate of 
incidence of the Blazhko effect in RRO stars in NGC 698 1 . We state 
that this is a lower limit because we are convinced that with more 
observations during our time-base and/or with an extension to our 
time-base, it is possible to detect more small amplitude Blazhko 
effects. 
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We also investigated secular period changes for each RR Lyrae 
star for which the light curve does not phase as well as expected at 
the best-fit period given our photometric precision, which included 
investigating the Blazhko effect variables discussed above. We used 
a modified version of the string-length method to search for the best 
fit of the following linear model for secular period change: 



m 



t - E 

T{t) 



t - E 

'ml 



P(t) = Po + Pit - E) 



(4) 
(5) 



where (f){t) is the phase at time t, P{i) is the period at time t, Pq is 
the period at the epoch E, and /3 is the rate of period change. For 
each light curve, the parameter space was searched at fixed epoch 
E, whose value is arbitrary, for the best-fit values of Pq and /?, 
using a brute force search for the minimum string-length statistic 
of the light curve, and using a small range of periods around the 
previously determined best-fit period. 

For the Blazhko effect variables described above, it is clear 
that the secular period change model cannot explain the observed 
light curve variations. However, for the variable V29, which we 
originally suspected of exhibiting the Blazhko effect, we find that 
the light curve produces a much better phased light curve for 
Po = 0.597472 d at epoch E = 2453507.446 d and with /3 ^ 
— 1.38xlO~* d d~^, indicating that the star pulsation frequency 
is slowly increasing over time. The derived value f or fi is rather 
large compared to typical values for RRO stars (e.g. ljurcsik et al.l 
l200lh although it is certainly plausible. 

In Figure [7] we show a portion of the V — r CMD for 
NGC 6981 that focusses on the horizontal branch and includes part 
of the red giant branch. As expected, the RRO stars cluster towards 
the red end of the instability strip, and the RRl stars cluster towards 
the blue end. There is, however, some scatter in the positions of the 
RR Lyrae stars, which is due to certain systematic errors in the 
measured mean magnitudes. Of the RR Lyrae stars with blended 
PSFs, V18 and V36 have clearly over-estimated mean magnitudes 
in V and r, which leads to a systematic error in their position on 
the CMD. The two other such stars, V24 and V51, do not suffer 
from this problem, presumably because the DanDIA software was 
able to perform successful deblending. Other outliers in the CMD 
may be explained by systematic errors associated with calculating 
the mean magnitude for an incomplete phased light curve (e.g. V2, 
V13, V25). 



3.6 Suspected RR Lyrae Stars 

Each of the variables V44, V45 and V53 lies close to a saturated 
star in the V and R reference images, and their light curves suf- 
fer from missing epochs caused by the saturated pixels (see Sec- 
tion |23](. In fact, for V44 and V45, the light curves only contain 
photometry from the iright of 20041005 (the night of best-seeing), 
and V44 is further missing photometry in the R filter. For V53, the 
light curve only contains photometry from the nights of 20041004 
and 20041005. It is this lack of data that has prohibited us from de- 
termining periods for these variables (see end of Section [373t . The 
light curves for V44, V45 and V53 are presented in Figure[8] 

The variable V52 has a very poor phase coverage with our ob- 
servations which has also prohibited us from determining a period. 
We present the light curve for this variable in Figure |9] 

From the position of V45, V52 and V53 in the instability strip 
of the CMD in Figure |7] it is obvious that these variables are RR 
Lyrae stars. We also believe that V45 is an RRO type variable be- 
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Figure 8. The V (black circles) and r (red circles) light curves for the vari- 
ables V44 (top panel), V45 (upper-middle panel) and V53 (lower two pan- 
els). The r light curves have been offset in magnitude such that the mean r 
magnitude matches the mean V magnitude for each star. Period determina- 
tion has not been possible due to the lack of data. 
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Figure 9. The V (black circles) and r (red circles) light curve for the vari- 
able V52. The r light curve has been offset in magnitude such that the mean 
r magnitude matches the mean V magnitude. Period determination has not 
been possible due to the veiy poor phase coverage of our data. 



cause the slow decline of the light curve over ~0.2 d is not con- 
sistent with a sinusoidal light curve of the period expected for a 
typical RRI type variable. We also believe that V44 is an RRO type 
variable because it has a mean V magnitude of ~ 17.04 mag, typ- 
ical of the RR Lyrae stars in this cluster, and it exhibits a decline 
in brightness similar to V45 over the same time window (see Fig- 
urelDl. 

3.7 SX Phoenicis Stars 

Recently, a large number of SX Phe type variables have been dis- 
covered in globular clusters in the blue straggler region of the 
CMD (e.g. R odrigue z & L6pez-Gonziie3 \20od . Ijeon et al.l l2003l: 
iDekanv & Kovacsl 120091) . In our search for new variables in 
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Figure 10. The V (black circles) and r (red circles) light curves for the SX Phe type variables: (a) V54, (b) V55, and (c) V56. The r light curves have been 
offset in magnitude such that the mean r magnitude matches the mean V magnitude for each star. Mean photometric uncertainties per data point in V and r 
are 15 and 15 mmag for V54, 31 and 35 mmag for V55, and 40 and 51 mmag for V56, respectively. 



NGC 6981, we have discovered three SX Phe type variables (V54- 
V56) and their light curves are presented in Figure[lO] We analysed 
each V and r light cu rve separately with the program PERIOD 04 
jLenz & Bregej2()05 ) in order to derive the dominant pulsation pe- 
riod and to search for any other pulsational frequencies. We ac- 
cepted any frequencies with oscillation amplitudes in V that ex- 
ceed the mean photometric uncertainties per data point in the cor- 



responding light curve (15 mmag for V54, 31 mmag for V55, and 
40 mmag for V56). The results of the frequency analysis for the V 
filter light curves are reported in Table |5] For V54, we detect two 
significant frequencies with frequency ratio /1//2 ~0.75, which 
therefore correspond to the fundamental (F) and first overtone (IH) 
radial oscillation modes. For V55 and V56, we did not detect any 
significant frequencies apart from the dominant frequency, and 
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Table 5. Detected pulsation frequencies for the SX Phe variables discovered 
in NGC 698 1 . For each SX Phe variable, we list the mean V magnitude Aq 
(column 2), the detected frequencies (columns 3 & 4), the corresponding 
(full) amplitude Ay in the V filter (column 5), and the identified mode 
of oscillation at each frequency (column 6). The numbers in parentheses 
indicate the uncertainty on the last decimal place. 



Variable 




Label 


Frequency 


Av 


Mode 


Star ID 


iV mag) 




(eye d^^) 


(mag) 




V54 


18.191(1) 


/i 


13.8915 


0.246 


F 






/2 


18.4544 


0.136 


IH 


V55 


19.392(3) 


/l 


21.2618 


0.174 




V56 


19.613(3) 


/l 


24.7099 


0.087 





therefore we are unable to identify the mode of oscillation. Note 
that the periods listed in Table [2] for these variables are the periods 
corresponding to the frequency of the largest amplitude oscillation. 

The positions of the variables V54-V56 in the calibrated V — I 
CMD are shown in Figure [12] as solid red circles. The / images 
were observed at one epoch only, and hence we have calculated the 
V and / mean magnitudes for each star from the three V images 
taken closest in time to the three / images, and from the three / 
images themselves, respectively. The V and / mean magnitudes 
differ in epoch by ~18 min, which is not important for the con- 
stant stars, but yields an approximate instantaneous colour for the 
variable stars. The fact that the V and / mean magnitudes are in- 
stantaneous rather than phase-averaged explains the spread of the 
RR Lyrae stars (solid and open blue circles) along the horizontal 
branch in Figure [12] and the fact that the SX Phe variable V54 lies 
close to, but not within, the dashed b ox delimiting the blue strag- 
gler region, which has been taken from lHarrii l ll993h for the cluster 
NGC 6366 and adapted to the distance and reddening of NGC 6981 
(see Section[53). The other two variables V55 and V56 lie within 
the blue straggler region. This evidence, along with the fact that the 
periods of the largest amplitude oscillation are in the range 0.03- 
0.08 d, confirms our SX Phe type classification for these variables. 

More observations of these variables are required to confirm 
the pulsational frequencies that we have detected, and to further 
characterise these stars. 



4 PHYSICAL PARAMETERS OF THE RR LYRAE STARS 
4.1 Fourier Light Curve Decomposition 

The Fourier light curve decomposition technique is based on fitting 
the following function to a light curve: 

m(t)^ Ao-^Y^^Ak cos \^—{t- E) + 4)k J (6) 

where m{t) is the magnitude of the object at time t, P is the period, 
E is the epoch, Ak and ipk are the fcth Fourier coefficients, and TV 
is the degree (or number of harmonics) of the Fourier decomposi- 
tion model. The Fourier parameters (relative amplitudes and phase 
differences), which are independent of E, are subsequently defined 
by the relations: 

R^j = A^/Aj (7) 

4'r] = J0» - i<t>] (8) 

for 1 < i,j < N. 
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Figure 11. Plot of <j>3i versus P for the RRO stars listed in Table [6] Open 
circles with labels represent the "outlier" RRO stars that were excluded from 
the fit of the cluster metallicity in Equations |9]to Equations 1 1 1 K continuous 
line). 



ISimon&Leel ( ll98ll) demonstrated that the light curve Fourier 
parameters could provide useful information about pulsating vari- 
ables, and currently, for Cepheids and RR Lyrae stars, the pa- 
rameters are routinely used to provide information about pulsation 
modes and resonances, and estimates of physical parameters such 
as metallicity, luminosity and effective temperature. With this mo- 
tivation, we have performed a Fourier decomposition of the V filter 
light curves of the RR Lyrae stars in NGC 6981 for which we have 
data and a reliable period estimate (i.e. 28 RRO and 5 RRl type 
variables). We note that the number of harmonics required for a 
good fit of a given light curve depends as much on the quality of 
the data, such as the amount of scatter and phase coverage, as on 
the light curve shape. Therefore we have chosen N for each light 
curve based on the highest degree harmonic for which the ampli- 
tude is significant. However, for some variables, the corresponding 
light curve does not yield a well constrained Fourier decomposi- 
tion due to a lack of full phase coverage, and we have abandoned 
the fits for the RRO type variables V2, V3, VI 1, V13, V21, V29 
and V32 due to this problem. For the RRO variables exhibiting the 
Blazhko effect, we used the same subsets of self-consistent data as 
those used for the period determination. 

In Table [6] we list the Fourier coefficients Ak for k = 
0, 1, 2, 3, 4, and the Fourier parameters (j!>2i, 031 and 041, for the 
21 RRO and 5 RRl type variables for which the Fourier decompo- 
sition fit was successful, along with the number of harmonics N 
used to fit each light curve. 

In previous papers (lArellano Ferro et al ] l2008al: 
lArellano Ferro et"aL I [2010) , we have discussed in detail the 
relations and corresponding zero points that may be used to 
calculate the metallicities [Fe/H], absolute magnitudes Mv, 
effective temperatures Tcff, masses M and radii R of the RR Lyrae 
stars in globular clusters from the fit parameters derived from the 
Fourier decomposition of the V light curves. Currently there is 
some debate about the reliability of these relations with respect 
to the calculation of e ffective temperatures, masses and radii (e.g. 
ICacciari. Corwin & C arnev 2005; Catelan 2004, see his Section 4). 
Therefore, in the following sections, we refrain from calculating 
masses and radii, and for the other physical parameters we simply 
proceed to derive their values for the RR Lyrae stars in NGC 698 1 
based on the published relations with little further discussion. 
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Table 6. Fourier coefficients A^^ for fc = 0, 1, 2, 3, 4, and the Fourier parameters 021, 4>3l and 04i, for the 21 RRO and 5 RRl type variables for which 
the Fourier decomposition fit was successful. The numbers in parentheses indicate the uncertainty on the last decimal place. Also listed are the number of 
harmonics N used to fit the light curve of each variable, and the deviation parameter Dm (see Section l4!2l . 



Variable Aq Ai A2 A-j, A4 <f>2i (f>3i 4>4.i N D, 

Star ID {V mag) {V mag) {V mag) {V mag) (V mag) 

RRO stars 



VI 


16.926(3) 


0.220(4) 


0.101(3) 


0.068(2) 


0.033(3) 


3.924(46) 


8.482(72) 


6.748(99) 


6 


2.7 


V4 


16.981(4) 


0.311(6) 


0.147(6) 


0.109(5) 


0.081(5) 


3.963(44) 


8.347(68) 


6.357(88) 


7 


2.1 


VS" 


16.969(32) 


0.355(41) 


0.167(46) 


0.150(35) 


0.098(23) 


4.480(359) 


8.211(486) 


7.022(667) 


7 


13.0 


Yl" 


17.036(13) 


0.334(21) 


0.204(17) 


0.152(13) 


0.063(13) 


3.762(110) 


7.455(153) 


5.602(215) 


6 


9.0 


V8 


16.966(6) 


0.291(8) 


0.141(11) 


0.101(10) 


0.063(6) 


3.945(73) 


8.254(109) 


6.273(181) 


8 


2.2 


V9 


16.941(4) 


0.236(5) 


0.102(6) 


0.072(5) 


0.037(4) 


3.958(67) 


8.417(99) 


6.670(147) 


5 


1.4 


VIO' 


16.971(10) 


0.302(17) 


0.143(9) 


0.115(8) 


0.075(6) 


3.836(72) 


8.018(81) 


6.004(136) 


8 


1.3 


VM* 


16.901(4) 


0.223(7) 


0.085(6) 


0.048(5) 


0.020(5) 


3.959(54) 


8.426(114) 


7.065(232) 


8 


2.6 


VIS*" 


16.987(11) 


0.314(15) 


0.169(15) 


0.105(14) 


0.074(11) 


4.105(134) 


8.316(198) 


6.234(276) 


10 


2.6 


V16 


16.905(3) 


0.284(4) 


0.134(5) 


0.098(4) 


0.060(4) 


3.933(41) 


8.249(63) 


6.349(81) 


7 


0.5 


V17 


16.968(5) 


0.295(3) 


0.147(6) 


0.103(7) 


0.062(6) 


4.004(84) 


8.305(117) 


6.392(170) 


10 


0.7 


VIS'^ 


16.376(58) 


0.296(97) 


0.115(60) 


0.056(58) 


0.031(39) 


3.979(731) 


8.086(900) 


6.278(1101) 


6 


2.6 


V20 


16.934(2) 


0.240(4) 


0.113(4) 


0.081(3) 


0.046(2) 


4.113(36) 


8.644(60) 


6.881(102) 


6 


1.4 


vas"^ 


16.948(2) 


0.241(3) 


0.142(4) 


0.109(4) 


0.051(3) 


4.177(34) 


8.611(47) 


6.678(82) 


6 


2.9 


V28<* 


16.968(4) 


0.266(5) 


0.133(5) 


0.071(5) 


0.031(4) 


4.108(35) 


8.633(74) 


6.960(169) 


8 


1.9 




16.958(20) 


0.303(19) 


0.181(23) 


0.065(23) 


0.063(16) 


3.509(265) 


7.590(411) 


5.773(464) 


5 


8.5 




16.503(5) 


0.214(6) 


0.124(6) 


0.062(5) 


0.040(5) 


3.746(76) 


8.229(140) 


6.266(179) 


4 


3.2 


V47 


16.926(1) 


0.112(2) 


0.030(1) 


0.011(2) 


0.001(2) 


4.286(71) 


8.669(126) 


8.513(1108) 


4 


2.0 


V48" 


16.894(4) 


0.176(5) 


0.074(5) 


0.038(3) 


0.018(3) 


3.955(92) 


8.591(168) 


6.935(241) 


5 


4.4 


V49* 


16.937(3) 


0.277(5) 


0.144(4) 


0.106(4) 


0.068(4) 


3.793(50) 


8.100(66) 


6.000(96) 


5 


2.7 


V50 


16.798(5) 


0.358(8) 


0.178(7) 


0.129(5) 


0.087(7) 


3.882(47) 


7.994(70) 


5.842(86) 


8 


2.4 



RRl stars 



V12 


16.856(2) 


0.251(2) 


0.047(3) 


0.020(2) 


0.018(2) 


4.667(46) 


2.594(122) 


1.044(131) 


7 — 


V24'= 


16.861(1) 


0.206(2) 


0.026(2) 


0.013(2) 


0.006(2) 


4.598(80) 


3.053(147) 


2.714(342) 


10 — 


V25 


16.882(3) 


0.204(4) 


0.030(5) 


0.021(3) 


0.009(3) 


5.359(83) 


3.726(205) 


2.427(454) 


5 — 


V43 


16.939(4) 


0.278(5) 


0.056(5) 


0.014(5) 


0.012(5) 


4.628(98) 


2.406(355) 


0.702(414) 


4 — 


V46 


16.806(1) 


0.149(2) 


0.009(2) 


0.005(2) 


0.003(2) 


4.498(196) 


3.519(367) 


0.284(695) 


4 — 



"These RRO variables have Dm > 3. ''These RRO variables exhibit a small amplitude Blazhko effect. '^Closely blended with stars of similar or greater 
brightness. ''These RRO variables exhibit a large amplitude Blazhko effect. 



4.2 Metallicity 

ljurcsik & Kovic3 ( Il996l) derived the following empirical 
[Fe/H]— </)3j' — P relation for RRO type variables based on 
a data set of V light curves and independent spectroscopic 
metallicity estimates for 81 field RRO variables: 

[Fe/H]j = -5.038 - 5.394P + 1.3450^"' (9) 

where (^3^ is a phase difference from the Fourier light curve de- 
composition using a sine series, as opposed to a cosine series, and 
P is the period (d). The value of (jif^ may be obtained from the 
equivalent phase difference 4>ij in a cosine series via: 

0!;' =<^.. (10) 

The metallicity [ Fe/H]j may be transformed to the ZW 
( IZinn& West! 1984) metallicity sc ale [Fe/H]zw by using the inverse 
of the relation from |jurcsikl ( ll995h : 

[Fe/H]j = 1.431[Fe/H]zw + 0.88 (11) 

The empirical relation defined in Equation |9] is only appli- 
cable to light curves of RRO variables that are "similar" to the 
light curves that were used t o derive Equation |9] For this pur- 
pose, ljurcsik & Kovacsl i 19961) defined a deviation parameter Dm 
describing the deviation of a light curve from the calibrating light 



curves, based on the Fourier parameter interrelations, and a com- 
patibility condition that Dm < 3. We list the values of Dm for our 
light curves in column 1 1 of Table |6] 

In column 2 of Table[7] we report the metallicity estimates for 
each RRO varia ble along with their unce rtainties as obtained via 
Equation 4 from ljurcsik & Kovacsl ( 1 19961) . The metallicities show 
a large spread of values ranging over ~1 dex due to a number 
of outliers. In order to highlight which RRO stars produce outlier 
metallicity values, we have plotted the parameters <^3i versus P 
in Figure [TT] for the RRO stars in Table [6] assuming that the un- 
certainty on P is negligible. We find that V7, V20, V23 and V28 
have clearly anomalous (j>-i,\ values given the associated uncertainty. 
However, this is somewhat to be expected since V23 and V28 ex- 
hibit large amplitude Blazhko effects, and V7 has a non-compatible 
light curve with the relation in Equation [9] (Dm > 3). For V20, we 
can find no reason why the light curve may yield an anomalous <j)zi 
value. 

We have decided to ignore the <^3i values for the RRO stars 
with non-compatible light curves (V5, V7, V31, V36, V48), for 
those exhibiting large amplitude Blazhko effects (V23, V28, V31), 
and for V20. These "outlier" variables are plotted in Figure [TT] as 
open circles. We note that it is not necessary to ignore the <^3i val- 
ues for the blended RRO stars since the 03 1 parameter is a light 
curve shape parameter that is unaffected by systematic errors in 
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Table 7. Physical parameters of the RRO variables. The numbers in paren- 
theses indicate the uncertainty on the last decimal place. 



Variable 


[Fe/H]zw 


Mv 


log(L/L0) 


Teff 


Star ID 




(mag) 




(K) 


VI 


— 1.45(10) 


0.601(5) 


1.669(2) 


6334(301 




— 1.33(9) 


0.623(8) 


1.660(3) 


6493(29) 




— 1.28(46) 


675r551 


1 639(221 


6483(2041 




—2 06f 16^ 


0.673(27) 


1.640(11) 


6344(641 


V8 


— 1 47tI2^ 


616(121 


1 663(51 


6433(471 


V9 


— 1.45(11) 


608(71 


I 666(31 


6358(41) 


VIO 


— 1.66(10) 


0.630(21) 


1.657(8) 


6408(36) 


V14 


-1.46(13) 


0.598(9) 


1.670(4) 


6308(50) 


V15 


-1.31(20) 


0.634(21) 


1.656(8) 


6519(84) 


V16 


-1.50(9) 


0.612(6) 


1.664(2) 


6408(26) 


V17 


-1.44(13) 


0.606(7) 


1.667(3) 


6430(49) 


VIS' 


-1.51(85) 


0.622(5) 


1.660(49) 


6441(372) 


V20" 


-1.21(9) 


0.622(5) 


1.661(2) 


6425(26) 




-1.20(8) 


0.657(5) 


1.646(2) 


6454(21) 


Yl?.°-'^ 


-1.11(10) 


0.622(7) 


1.660(3) 


6476(34) 




-2.00(39) 


0.611(29) 


1.665(12) 


6324(165) 


V36"''' 


-1.55(15) 


0.653(8) 


1.648(3) 


6366(57) 


V47 


-1.39(14) 


0.641(3) 


1.653(1) 


6119(117) 


V48'' 


-1.43(17) 


0.601(6) 


1.669(3) 


6294(69) 


V49 


-1.65(9) 


0.623(7) 


1.660(3) 


6388(28) 


V50 


-1.42(9) 


0.684(10) 


1.635(4) 


6569(30) 


Weighted Mean: 


-1.48(3) 


0.623(2) 


1.660(1) 


6418(10) 



"These variables are not included in the calculation of the mean 
metallicity (see Section |4!2l . ''These variables are not included in the 
calculation of the mean absolute magnitude, log-luminosity or effective 
temperature (see Sections |4.3l & |4~4l . '^V4 and V23 have spectroscopic 
metallicity measurements of — 1.61 and —1.28, respectively, on the ZW 
scale (see Section l4!2t . 



Table 8. Physical parameters of the RRl variables. The numbers in paren- 
theses indicate the uncertainty on the last decimal place. 



Variable 


[Fe/H]zw 


Mv 


log(L/LQ) 




Star ID 




(mag) 




(K) 


V12 


-1.59(21) 


0.568(9) 


1.651(4) 


7373(35) 


V24 


-1.77(27) 


0.587(10) 


1.643(4) 


7280(32) 


V25 


-1.70(42) 


0.515(14) 


1.672(6) 


7262(229) 


V43 


-1.63(59) 


0.601(23) 


1.638(9) 


7372(74) 


V46'" 


-1.04(66) 


0.643(12) 


1.621(5) 


7466(54) 


Weighted Mean: 


-1.66(15) 


0.568(6) 


1.651(2) 


7327(22) 



"This variable is anomalous in the Bailey diagram (see Figure|6) and it has 
an anomalous </)3i value. Therefore, we have not included this variable in 
the calculation of the mean metallicity, absolute magnitude, log-luminosity 
or effective temperature (see Sections [4.2ll4.3l & l4.4t . 



the light curve amplitude due to blending (see Section [23}. We 
used the remaining RRO stars to fit a metallicity [Fe/H]zw un- 
der the assumptions that the relation in Equation |9] is correct and 
that NGC 6981 has no intrinsic spread in metallicity. We obtain 
[Fe/H]zw ~ — 1.48±0.03 and plot the corresponding relation as 
a straight line in Figure [TT] Clearly the ljurcsik & Kovacsl h996l) 
empirical [Fe/H]— (^3/ — P relation holds for the RRO variables in 
NGC 6981. It is also satisfying to note that only the stars V7, V20, 
V23 and V28 have 031 values with residuals greater than 3(t. 

For the RRl variable s, we used the empirical [Fe/H]— 03 i — 
P relation determined by [Morgan. Wahl & Wieckhor^ j2007h in 



order to estimate metallicities: 

[Fe/H]zw =52.466P^ - 30.075P -f 0.131(^3i 
+ 0.982</)3i - 4. 198(^31 P + 2.424 

(12) 

This relation provides metallicities on the ZW scale with a standard 
deviation in the residuals of 0.145 dex for the 106 cluster stars in 
the calibrating sample. The [Fe/H]zw values for the RRl variables 
are reported in column 2 of Table [8] along with their uncertainties 
calculated from the propagation of the (^31 uncertainties, but ignor- 
ing any uncertainty in the period. 

We note that for the RRl variable V46, the value of the 031 
parameter is anomalously large given the short period, and this has 
resulted in the derivation of an anomalous metallicity of —1.04 
for this star. In the Bailey diagram in Figure |6] we have already 
noted that V46 also has an anomalously small amplitude for its pe- 
riod, and therefore we have decided to exclude this variable from 
the calculations of the mean physical parameters of the RRl vari- 
ables in Table [8] Consequently, we calculate a mean metallicity of 
-1.66±0.15 for the RRl stars. 

Finally, we note that the metallicities of only two RR Lyrae 
stars in NGC 6981 have previously been measured in the literature. 
ISmith&Perkin3 ( ll982[) measured the spectroscopic metallicity in- 
dex AS* for the RRO type variables V4 and V23 as 7.6 and 5.5, re- 
spectively, which yield metallicities [Fe/H]zw ~ —1.61 and— 1.28, 
respectively , when converted to the ZW scale via the following re- 
lation from lSuntzeff. Kinman & KraftI ( Il99lh : 

[Fe/H]zw = -0.158AS - 0.408 (13) 

We include these metallicity measurements as a footnote to Table|7] 
for the sake of completeness. 



4.3 Absolute Magnitude And Luminosity 

The following empirical relation yielding the absolute magnitude 
Mv of the RRO variables w as determined from V lig ht curve data 
for 383 such cluster stars bv lKovacs & Walkeil ( [200i1) : 

Mv = -1.876 log P- 1.158^1 + O.82IA3 + A'o (14) 

where Kq is a constant to be de t ermin ed. We adopt Ko — 
0.41 mag, as in [Arellano Ferro et al. in order to be consis- 

tent with a true distan ce modulus for the LMC of fio =18.5 mag 
jFreedman etai]|200lh . We list the derived values of Mv for the 
RRO variables in co lumn 3 of Table |7] 

lKovacsl ( ll998l) derived a similar empirical relation for Mv for 
RRl variables from V light curve data for 93 cluster stars: 

Mv = -0.961P - 0.0440^1' - 4.447A4 + Ki (15) 

where Ki is a constant to be determined. A value of Ki — 
1.261 mag was calculated by Kovacs (1998) to place Equa- 
tion 1151 on the Baade-Wesselin k luminosity scale. However, 
ICacciari. Corwin & CarnevI ( l2005h have suggested that this value 
should be decreased to Ki = 1.061 mag to make Equation|15|con- 
sistent with a LMC true distance modulus of ^0 ~ 18.5 mag. We 
have therefore adopted Ki = 1.061 mag, and we list the derived 
values of A'lv for the RRl variables in column 3 of Table[8] 

The empirical relations in Equations[T4]&[T5]depend on a sub- 
set of the amplitude Fourier coefficients. Therefore, any variable 
which suffers from a systematic over-estimate of the reference flux 
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due to blending will yield Fourier amplitudes that are systemati- 
cally too small, which in turn leads to a systematic error in the cal- 
culated value for Mv- Of the RR Lyrae variables with Fourier de- 
composed light curves, VI 8, V24 and V36 are blended with stars of 
similar or greater brightness, although only V18 and V36 have no- 
ticeably brighter mean magnitudes Ao than the rest of the RR Lyrae 
stars (see Table[6l(. Consequently, we ignore the Mv values for V18 
and V36 in calculating the mean absolute magnitude in Table [7] 
We also ignore the My values for V23, V28 and V31 because the 
strong Blazhko effect in these variables is likely to render the light 
curve amplitude measurements as unreliable. Finally, we note that 
the compatibility condition described in Section |4~2l onlv applies to 
the empirical relation in Equation |9] Hence we calculate mean ab- 
solute magnitudes of Mv ^ 0.623±0.002 and 0.568±0.006 mag 
for the RRO and RRl variables, respectively. 

We converted the absolute magnitudes listed in Tables |7]& [8] 
to log-luminosities log(L/L0) using: 



log (L/Le 



-0.4 {Mv + Bc(reff) - Mboi,©) 



(16) 



where we adopt a bolometric solar absolute magnitude of Mboi,© = 
4.75 mag, and a bolometric correction -Bc(reff) as a function of 
effec tive temperature from th e calibrations for metal-poor stars 
from lMontegriffo et alj ( Il998h . Anticipating our results from Sec- 
tion 14.41 we use the mean effective temperatures of the RRO 
and RRl variables presented in Tables |7] & [8] to determine bolo- 
metric correction s of -0 .023 and 0.055 mag, respectively, from 
iMontegriffo et alj ( Il998h , where we have reduced the mean effec- 
tive temperature of the RRl stars by 250 K in order to bring the 
mean temperature in line with temperatures derived from observed 
V — I colours (i Arellano Ferro et al.201Cl) . 



4.4 Effective Temperature 

The effective temperature Tcff of RRO and RRl variables may also 
be estimated via the Fourier parameters of the V ligh t curv es. For 
the RRO variables, we used the relations from JurcsiH h998l) based 
on a comparison of empirical data to the lKuruczl ^1993h model at- 
mospheres: 

(V - K)o = 1.585 + 1.257P - 0.273Ai - 0.2344"' 

+ 0M2(t>i]> (17) 
log Teff =3.9291 - 0.1112(1/ - K)o - 0.0032[Fe/H]j (18) 



in column 5 of Table [S] and calculate a mean effective temperature 
for the RRl variables of Teff ^ 7327±22 K. 

However, we urge caution in the use of the effective tem- 
peratures derived from the F ourier decomposit ion fit para meters 
since it has been shown by ICacciari. Corwin & Carnev l l2005l) 
that the temperatures calculated via Equations [TS] & [T9] do not 
match the co lour-temperature relations pre dicted by the tempera- 
ture scales of Sekiguchi & Fukugital ([2000) or by the evolutionary 
models of Cast elli (1999). In our previous works on NGC 5466 
dArellano Ferro e t al. 2008a) and NGC 5053 ^Arellano Ferro et all 
|2010|) . we also found similar systematic errors in the effec- 
tive temperatures derived from the light curve Fourier parame- 
ters for the RRl stars. Ho wever, for the RRO stars, we found 
dArellano Ferro et^ l201(]h that the effective temperatures de- 
rived from the light curve Fourier parameters agree with the 
colour-temperatures predic ted by t he horizontal branch models of 
VandenBerg, Bergbusch & Dowlerl |2006). Despite these caveats, 
we still provide the effective temperatures for the RR Lyrae stars 
in NGC 6981 in order to enable comparisons with similar work for 
other globular clusters. 



5 CLUSTER PHYSICAL PARAMETERS 
5.1 Oosterhoff Classification 

We calculate the mean periods for the RRO and RRl stars in 
NGC 6981 as (Prm) =0.563 d and (Prr,) =0.308 d, respec- 
tively. We have used the periods derived in this work for all RRO 
and RRl variables except V 51 due to its unreliable period, and 
we included the periods from lDickens & FlinnI d 19721) for the RRO 
variables V27 and V35 which do not lie in our field of view. We 
also calculate the ratio of the number of RRl to RR Lyrae type 
variables as tirri/ (wrro + riRRi) «0.14 using 36 RR Lyrae stars 
with robust classifications (i.e. excluding V44, V45, V52 and V53). 
These quantities are clearly consistent with the cl assific ation of 
NGC 6981 as an Oosterhoff type I cluster ( Smith 1995^, which 
agrees with previous cla ssifications ( ,van Albada & Bakeij 
ICastellani & Ouartal [19871) based on the jPickens & FlinnI 
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survey of RR Lyraes, but employs an updated and more complete 
catalogue of RR Lyrae stars for a more robust result. 



where [Fe/H]j is the metallicity calculated via Equation [9| We list 
the derived values of Teff for the RRO variables in column 5 of Ta- 
ble |7| 

EauationI 1 7 [depends on the amplitude Fourier coefficient Ai, 
and therefore, for the same reasons given in Section [43] we ig- 
nore the variables VI 8, V23, V28, V31 and V36 in the calcu- 
lation of the mean effective temperature for the RRO variables 
Teif «6418±10 K. 

For the RRl stars, one may use the following relation derived 
from the Fourier decomposition of theoretical light curves gener- 
ated from hydro-dynamical pulsation models dSimon & Clemend 
Il993h to estimate the effective temperature: 

log Teff = 3.7746 - 0.1452 log P + 0.0056(^3i (19) 

Note that the temperatures calculated from this equation are on 
a different absolute scale than the temperatures calculated from 
Equationll8l We list the derived values of Tcff for the RRl variables 



5.2 Cluster Metallicity 

According to lSandag^ d2006l ). the mean log-period of the field RRO 
type variables is related to metallicity on the ZW scale as follows: 

(log Prro) = -0.416 - 0.098[Fe/H]zw (20) 

[element et al. [ d200lh also previously observed that such a corre- 
lation exists for the Oosterhoff type I clusters, but that it breaks 
down for the Oosterhoff type II clusters. However, since NGC 6981 
is of Oosterhoff type I, we may employ Equation [20| to derive 
[Fe/H]zw ~ —1.68 based on (log Prro) ~ —0.251 d calculated 
for the same RRO stars as in Section [STl and assuming that the 
relation also applies to cluster RR Lyraes. 

The shortest and longest RRO periods depend on the loca- 
tion of the blue and red edg es, respectively, of the RRO instability 
strip, whichlSandagd 11993') showed also depend on the metallicity. 
ISandagd |T993) quotes the following relations between metallicity 
and the shortest P^in and longest Pmax periods for globular cluster 
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Table 9. Metallicity estimates for NGC 698 1 on the ZW scale as found from an extensive literature search. 



lFe/Hl7w 


Reference 


Method 


-1.48±0.03 


This work 


Fourier light curve decomposition of the RR Lyrae stars 


-1.46 


^araston et al. (2003J 


Lick indices 


-1.41 ±0.30 


Thomas, Maraston & Bender (2003) 


Lick indices 


-1.26±0.07 


Geisler. Claria & Minniti (1997) 


Washington photometry 


-1.50±0.05 


Rutledee. Hesser & Stetson (1997) 


Ca 11 triplet 


-1.4±0.1 


Brocato etal. (1996) 


CMD analysis 


-1.40 


Harris (1996) 


Globular cluster catalogue" 


-1.72±0.19 


Brodie&Huchra(1990) 


Absorption-line indices 


-1.5 


Rodgers & Harding (1990) 


Calcium abundance 


-1.54±0.09 


Zinn & West (1984) 


Qag spectral index 


-1.58±0.12 


Zinn (1980) 


Qag spectral index 



"^The catalogue version used is the updated 2003 version available at http://www.physics.mcmaster.ca/Globularhtn3] 



RRO stars (his Figure 10): 

log Pmm = -0.526 - 0.117[Fe/H]zw (21) 
log Pmax = -0.248 - 0.070[Fe/H]zw (22) 

In NGC 6981, V2 and V27 are the RRO variables with the shortest 
(-Pm,„ = 0.465254 d) and longest known (Pn,ax = 0.673774 d) pe- 
riods, respectively, from which we derive metallicities of ~ —1.66 
and ~ —1.09 via Equati ons 1211 & [221 r espectivelv. The updated 
version of Eg uation 1 2 1 1 in I S andagd ( |2006[) also yields a metallicity 
of ~ —1.66 after taking into account the period shift of the short- 
period locus between field and cluster RRO variables. 

The large discrepancy between the two metallicity values we 
have derived from considering the location of the blue and red 
edges of the instability strip would be reduced if any of the vari- 
ables V44, V45, V51, V52 and V53 are found to be RRO stars with 
periods shorter than 0.465 d or longer than 0.674 d. In fact, Equa- 
tion[22]is the most sensitive to a change in the period boundaries, 
and all that is required to bring the metallicity estimates from the 
shortest and longest periods into agreement at [Fe/H]zw ~ —1.66 
is a ~ 10 per cent increase in the value of P^ax to ~0.74 d. We 
also note that any future periods determined for V44, V45, V51, 
V52 and V53 will influence the mean (log-)period of the RRO vari- 
ables in NGC 6981, which may also have a substantial effect on the 
metallicity estimated via Equation|20|due to the high sensitivity of 
this relation on the mean log-period. In fact, we remain rather scep- 
tical of the metallicity estimates for NGC 698 1 that we have derived 
via Equations |20||21l &l22lfro m the work of Sandage for just these 
reasons, and Figure 4 from Ca cciari. Corwin & Carney ( 2005) pro- 
vides a clear graphical illustration of the potential problems with 
this method that we have described. 

We have compiled estimates of [Fe/H] on the ZW scale given 
in the literature for NGC 698 1 and obtained using a variety of meth- 
ods, and we summarise the results in Table[9l which includes refer- 
ences and a very brief description of the methods employed. 

In Section 14.21 we calculated mean metallicities of 
-1.48±0.03 and -1.66±0.15 on the ZW scale for the RRO and 
RRl type variables, respectively. The smaller uncertainty on the 
result for the RRO stars is mainly a consequence of the inclusion 
of more stars in the calculation of the mean metallicity for the RRO 
variables (13 RRO stars as opposed to 4 RRl stars). It should be 
noted that the uncertainties we have quoted only represent the in- 
ternal error for each metallicity estimate, and do not include any 
systematic errors that may be inherent in the Fourier light curve de- 
composition method that we have employed. Given the fact that dif- 
ferent empirical relations have been used to estimate the metallici- 



ties of the RRO and RRl stars, there may be some systematic offset 
between the metallicity estimates of the two types of variable. How- 
ever, assuming that this is not the case, and performing a combined 
mean metallicity estimate for all 17 stars also yields a mean metal- 
licity of ~ -1.48±0.03. Hence we adopt [Fe/H]zw ~ -1.48±0.03 
as our estimate of the metallicity of NGC 6981 on the ZW scale 
via the Fourier decomposition of the light curves of the RR Lyrae 
variables, which lies in the middle of the distribution of the other 
metallicity estimates for the cluster from a large range of indepen- 
dent methods th at are listed in Table [9l 

Recently, Carretta et ^([200^ collected a homogeneous set 
of intermediate to high-resolution spectra of ~2000 red giant 
branch stars in 19 Galactic globular clusters and used them to de- 
fine an accurate and updated metallicity scale. They also provide 
the following transformation from the ZW scale to this new metal- 
licity scale: 

[Fe/H]uvES = -0.413 + 0.130[Fe/H]zw - 0.356[Fe/H]^w (23) 

where [Fe/H ]uvES represents the metallicity on the new 
ICarretta et al. I {2009) scale. Our estimate of the metal- 
licity of NGC 6981 transformed onto this new scale is 
[Fe/H]uvES ~ — 1.38±0.03. This val ue is consistent to wi thin the 
uncertainties with the metallicity that lCarretta et alj ^20091) derive 
for NGC 6981 of [Fe/H]uvES ~ -1.48±0.07, which is based on a 
recalibration of the (Jsg and W' spectral indices, and a weighted 
mean of the metallicity estimates available on different scales, but 
transformed to the new scale. 

5.3 Distance 

To determine the distance to NGC 6981, the RR Lyrae stars can 
serve as standard candles since we have estimated their mean 
absolute magnitudes in Section 14.31 as Mv ~ 0.623±0.002 and 
0.568±0.006 mag for the RRO and RRl variables, respectively. We 
calculate the average of the RRO and RRl star mean V magnitudes 
(Ao from Table|6l( as -16.927±0.001 and ~16.865±0.001 mag, 
respectively (excluding the variables V18, V23, V28, V31 and V36 
from the calculation for the RRO stars because of the likely influ- 
ence of blending and the Blazhko effect on the measured mean 
magnitudes for these stars, and excluding the variable V46 from the 
calculation for the RRl stars because of its anomalous status). This 
implies mean apparent distance moduli of ~ 16.303±0.002 and 
16.297±0.006 mag for the RRO and RRl variables, respectively. 

NGC 6981 is known to be a globular cluster with relatively 
low interstellar extinction. The many values in the literature for 
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Table 10. True distance moduli (column 1), and distance estimates (column 2), for NGC 6981 as found from an extensive literature search. 



/iO (mag) 


Distance (kpc) 


Reference 


Method 


16.117±0.047 


16.73±0.36 


This paper 


Fourier light curve decomposition of the RRO stars 


16.111±0.047 


16.68±0.36 


This paper 


Fourier light curve decomposition of the RRl stars 


16.42±0.01 


19.23±0.09 


Cassisi. De Santis & Piersimoni (2001) 


Magnitude of the zero-age horizontal branch 


16.52 


20.14 


ZoccaU & Piotto (20001 


Globular cluster luminosity function 


16.08 


16.44 


Ferraro et al. (1999) 


Magnitude of the zero-age horizontal branch 




17.0 


Harris (1996) 


Globular cluster catalogue" 


16.30±0.05 


18.20±0.42 


Jimenez et al. (1996) 


CMD analysis 




16.4 


Zinn (1980) 


Q39 spectral index 


16.29 


18.1 


Searle & Zinn (1978) 


Spectral indices 


16.30 


18.2 


Woltier(1975) 


Mean magnitude of the RR Lyrae stars 


15.73 


14.0 


Dickens & Flinn (1972) 


Mean magnitude of the RR Lyrae stars 


17.0 


25 


Kron & Mavall (1960) 


Mean magnitude of the RR Lyrae stars 



"The catalogue version used is the updated 2003 version available at |http://www.physics.mcmaster.ca/Globular.htmI| 



the associ ated interstellar redd ening range from E{B — V) = 
0.00 mag (Searle & Zinn 1978) to E{B ~ V) =0.11 ±0.03 mag 
dRodgers & Harding. 1990) . with a median value of E{B — V) = 
0.06 mag. The two most recent determinations include E{B — 
V) = 0.06 mag from the analysis of 100 / xm emission from inXsr- 
stellar dust Jschlegel, Finkbeiner & Davy 1998). and E{B~V^ = 
0.070±0.015 mag from the comparison of predicted and measured 
colours as a function of i? light curve amplitude for the RRO stars 
jPe Santis & Cassisilll999l) . For the purpose of our analysis, we 
adopt a reddening of E{B -V)= 0.06±0.015 mag. 

Assuming Rv =3.1 for our Galaxy, we derive mean true dis- 
tance moduh of ^0 « 16.117±0.047 and 16.111±0.047 mag for 
the RRO and RRl variables, respectively, which translate to mean 
distances of ~16.73±0.36 kpc and ~16.68±0.36 kpc, respectively. 
Even though both of these distance estimates are consistent with a 
true distance modulus for the LMC of ^lq = 18.5 mag, they are 
calculated via two different empirical relations, each of which has 
its own systematic uncertainty, and therefore we cannot simply cal- 
culate an overall distance estimate from all RR Lyrae stars put to- 
gether. 

In Table [Tol we report our distance estimates for NGC 6981 
along with those we found in the literature search. We list both true 
distance moduli (where available) and distances (kpc) along with 
references and methods employed. The distances from the literature 
show considerable spread, although varying assumptions about the 
amount of reddening towards NGC 6981 certainly contribute sig- 
nificantly to the scatter. Our distance estimates would range from 
~15.5 kpc to ~18.2 kpc by simply changing our assumption about 
the reddening to values in the range 0.00 < E{B — V) <0. 11 
(defined by the minimum and maximum measured reddening val- 
ues in the literature). We observe that our distance estimates for 
NGC 6981 fall towards the shorter end of the distribution of dis- 
tance estimates from the literature, and it is satisfying to note that 
our mean distances for both types of RR Lyrae star are in excellent 
agreement. 



5.4 Age 

It is only in the last 20 years that measurements of the age of 
NGC 6981 have appeared in the literature, and they number far 
fewer then the number of measurements of the metallicity and dis- 
tance. We present the measured ages from the literature in Tablellll 
and we observe that the measurements have converged to a value 
of ^^13 Gyr. 



An exquis itely precise and deep C MD for NGC 698 1 has been 
constructed by lAnderson et al.l ([2008') through careful processing 
of a stack of images of NGC 698 1 obtained with the Wide-Field 
Channel of the Advanced Camera for Surveys (ACS) on board the 
Hubble Space Telescope. These images were observed as part of 
the ACS Survey of Globular Clusters which includes observations 
of the central ~3 x3 arcmin^ of 65 Galactic globular clusters. Sub- 
sequent analysis of all the cluster CMDs in the survey via main se- 
quence fitting and main-sequence turnoff (MSTO) measurements 
allowed |m arfn-Franch et al.l (I2OO9') to determine pre cise (2%-7% 
uncertainty) relative ages between the clusters. Dotter "et"ai] ( 1201 Cl) 
placed these relative ages on an absolute scale, thereby deriving an 
age of 12.75±0.75 Gyr for NGC 6981, which is the most precise 
age measurement to date. 

Our calibrated photometric data for NGC 6981 in the V and I 
bands does not reach as deep as the corresponding ACS data. Also, 
in our observations, most stars in the central region (3x3 arcmin^) 
of the globular cluster are blended with other cluster members, 
making accurate absolute photometry a very difficult task for the 
vast majority of cluster members, and introducing systematic pho- 
tometric errors due to poorly modelled PSFs and blended light con- 
tamination. The space-based observations enjoy a much better spa- 
tial resolution than our data due to the lack of blurring from an 
intervening atmosphere, and therefore the ACS photometry is of 
much better accuracy with considerably less systematic errors due 
to blending. Hence, any analysis of the CMD from our data for de- 
ter mining the c l uster age cannot compete with the similar analysis 
bv lDotter et ai] ( I2OI0I) of the ACS CMD. 

In the calibrated V — I CMD presented in Figure[T2l we over- 
plot a theoretical isochrone (solid black line) corresponding to the 
iDotter et al.l ([2010) age of 12.75 Gyr, our derived metallicity of 
[Fe/H]zw = —1.48, an a-element abundance of [a/Fe]= +0.3 
(a typical assumption for globular clusters), and an apparent dis- 
tance modulus of ^ = 16.303 mag (our result for the RRO stars). 
The isochrone has been linearly interpolated from the four near- 
est isochrones in age and metallicity as taken from the stella r 
evolutionary models of lVandenBerg. Bergbusch & DowleJ ( |2006|) . 
which are also plotted in Figure [T2| as solid coloured lines. The 
isochrones show that our CMD i s consistent with the age for 
NGC 6981 derived bv lDotter et al] ([TolO), although the fit to the 
upper part of the red-giant branch is too red by about ~0.1 mag. 
This mismatch can be slightly improved by adopting a more metal- 
poor isochrone but at the expense of a worsening fit near the main 
sequence turn-off. Note that we have not applied a.V ~ I colour 
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Age (Gyr) 


Reference 


Method 


12.75±0.75 


Dotter et al. (20101 


Isochrone fitting to the CMD 


9.5 


MeisMLejL& Weiss (2006) 


Isochrone fitting to the CMD 


12.6±1.0 


Rakos & Schombert (2005) 


Stromgren photometry 


14 


Zoccali & Piotto (2000) 


Globular cluster luminosity function 


12.7±2, 13.0±2, 13.5±2 


Jimenez et al. (1996) 


CMD analysis 


12.4±2.2 


Saraiedini & Kins (1989) 


CMD analysis 



14 5 L[Fe/H]=-1.48 ; age= 12.75Gy (Black) 




V I 

Figure 12. The V - / CMD for NGC 6981. Note that the magni- 
tude measurements approximate instantaneous magnitudes which explains 
the spread in the positions of the RR Lyrae stars along the horizon- 
tal branch (solid and open blue circles for the RRO and RRl variables, 
respectively) and the fact that one of the SX Phe variables (solid red 
circles) does not lie within the dashed box delimiting the blue strag- 
gl er region. Theoretical isochrones f rom the stellar evolutionary models 
of IVandenBerg. Bergbusch & DowleJ ^2006.) are plotted as solid coloured 
lines. The black solid isochrone has been linearly interpolated from the 
other isochrones to the iDotter et"^ ] i2oToh age of 12.75 Gyr and our de- 
rived metallicity of [Fe/H]zw = — 148. 

correction to the isochrones due to reddening since this correc- 
tion makes the isochrone fit considerably worse by shifting the 
isochrones by ~0.076 mag to the right (assuming E{B — V) = 
0.06 mag). This implies that the adopted reddening is too large 
and/or that the adopted theoretical isochrones are systematically 
too red at the cluster age and metallicity. 



6 SUMMARY AND CONCLUSIONS 

We employed the technique of difference image analysis to per- 
form precision differential photometry on a set of CCD time-series 
observations of the globular cluster NGC 6981 in the Johnson V , 
R and I filters. The difference imaging technique has allowed us to 
study in detail the variable star population of the cluster, including 
in the highly crowded central regions where traditional PSF fitting 
fails. Compared to previous photographic studies of this cluster, 
our photometry is deeper by ^^4 magnitudes, achieves a much bet- 
ter precision per data point (<20 mmag down to ~18.5 mag as 



opposed to ~50 mmag at ~17 mag), and covers a longer time-base 
of ~5 years with substantially more data points (~100 compared 
to ~20-60). 

Consequently, we have been able to perform a census of the 
variable star population in NGC 698 1 , and to clarify the status of al- 
ready known and suspected variables in the cluster We have shown 
that 20 of the 58 stars that are labelled as variables in the litera- 
ture are actually non- variable, and we present the details of all 43 
confirmed variables in NGC 6981 in Table|2] 

We have discovered 14 of the 43 confirmed variables in 
NGC 6981 from our data. The new detections consist of 11 RR 
Lyrae variables and 3 SX Phe variables, although we were only able 
to derive reliable periods for 6 of the newly discovered RR Lyrae 
stars. For the 27 of the 29 previously known RR Lyrae variables that 
lie in our field of view, we calculated a new set of ephemerides that 
improve substantially on their previous values in terms of accuracy. 

The current set of 43 variables in NGC 6981 consists of 40 
RR Lyrae stars and 3 SX Phe stars. Furthermore, the set of RR 
Lyrae stars is made up of 31 RRO type variables, 5 RRl type 
variables, and 4 RR Lyrae variables with an ambiguous subtype, 
although these variables are most likely to be of the RRO type 
(V44, V45, V52 and V53). Based on the mean periods for the RRO 
and RRl variables, and the ratio of the number of RRl to RRO 
type variables, we confirm the Oosterhoff type I classification for 
NGC 6981. We report the detection of a strong Blazhko effect in 5 
RRO variables (VI 1, V23, V28, V31 and V32), and a smaller am- 
plitude Blazhko effect in another 5 RRO variables (VIO, V14, V15, 
V36 and V49), which implies a ~34 per cent lower limit on the 
rate of incidence of the Blazhko effect in RRO stars in NGC 6981. 
We also find that the RRO star V29 exhibits a secular period change 
of /? « — 1.38x 10~* d d~^, indicating that the star pulsation fre- 
quency is slowly increasing over time. 

Our analysis of the light curves of the SX Phe variables has 
allowed us to detect two pulsation frequencies in the variable V54, 
corresponding to the fundamental and first overtone radial oscilla- 
tion modes. We did not detect any other frequencies apart from the 
dominant frequency for the other two SX Phe variables, and there- 
fore we are unable to identify their modes of oscillation. 

We provide calibrated V and / photometry in the Johnson- 
Ki'on-Cousins photometric system, derived from the analysis of a 
set of standard stars in the field of the cluster, and instrumental 
r photometry, for all of the 41 confirmed variables in NGC 6981 
in our field of view. This data is available in electronic form (see 
Supporting Information) and the format is shown in Table|4] In or- 
der to collect the information required to observe these objects into 
one place, we calculate celestial coordinates for each variable (see 
Table [3]( and provide detailed finding charts in Figure [3] Finding 
charts for the two RR Lyrae varia bles V27 and V 35 that are not in 
our field of view may be found in lDickensl ( [T972l) . 

We performed a Fourier decomposition of the light curves for 
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21 RRO and 5 RRl variables with reliable period estimates and suit- 
able phase coverage, and we report the corresponding Fourier coef- 
ficients in Table |6l The Fourier parameters have been used to esti- 
mate the metallicity, absolute magnitude, log-luminosity, and effec- 
tive temperature for each RR Lyrae star based on empirical, semi- 
empirical, and theoretical relations available in the literature. As- 
suming that the RR Lyrae stars in NGC 6981 are of the same com- 
position, distance and age, then appropriate averages of the derived 
properties of the RR Lyrae stars may be calculated and employed 
as estimates of these properties for the parent cluster. Applying this 
method, we derive a metallicity of [Fe/H]zw ~ — 1.48±0.03 on 
the ZW scale for NGC 6981 , and [Fe/H]uvES ~ -1.38±0.03 on 
the new ICarrettaetalJ (|2009i scale. Similarly, we derive mean true 
distance moduli of /xo « 16. 1 17±0.047 and 16. 1 1 1 ±0.047 mag for 
the RRO and RRl variables, respectively, and corresponding dis- 
tances to NGC 6981 of ~16.73±0.36 kpc and ~16.68±0.36 kpc, 

respectively. 

The age of NGC 6981 has been estimated by iDotter et al.l 

1 I2OICII) as 12.75±0.75 Gyr from analysis of Hubble Space Tele- 
scope ACS observations, and our calibrated V — I CMD cannot 
compete with their data in terms of accuracy and lack of systematic 
errors. Hence we simply demo nstrate that our CMD data is consis- 
tent w ith the isochrones of the IVandenBerg. Bergbusch & DowleJ 
( l2006h stellar evolutionary models when interpolated to the age and 
metallicity of NGC 6981, and if one chooses to ignore the colour 
effect of the reddening. It is likely that the reddening is smaller 
than we assumed in Section [53] or that there is a systematic error 
in the colour of the theoretical stellar evolutionary models that we 
adopted, or both. 
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Table 4. Time-series V, r and / photometry for all the con- 
firmed variables in our field of view, except V27 and V35 which lie 
outside of our field of view. 
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tent or functionality of any supporting materials supplied by the au- 
thors. Any queries (other than missing material) should be directed 
to the corresponding author for the article. 
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